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ABSTRACT 

Context. The Pierre Auger Collaboration has reported 27 Ultra-High Energy Cosmic Ray Events (UHECRs) with energies above 
56xl0 18 eV (56 EeV) and well determined arrival directions as of 2007 August 31. They find that the arrival directions of these 
UHECRs are not isotropic, but instead appear correlated with the positions of nearby Active Galactic Nuclei (AGNs) from the catalog 
of Veron-Cetty & Veron. 

Aims. Our aim was to determine the sources of these UHECRs by comparing their arrival directions with more complete and/or 
comprehensive astronomical source catalogs. 

Methods. We have cross correlated the arrival directions of the UHECRs with the positions of supernovae, radio supernovae, galaxies, 
active galaxies, radiogalaxies, and clusters of galaxies, all at distances within ~200 Mpc. 

Results. Four (eight) of the 27 UHECRs with energy greater than 56 EeV detected by the Pierre Auger Observatory have arrival 
directions within 1.5° (3.5°) of the extended (> 180 kpc) radio structures of nearby radiogalaxies or the single nearby BL Lac with 
extended radio structure. Conversely the radio structures of three (six) of all ten nearest extended radiogalaxies are within 1.5° (3.5°) 
of a UHECR; three of the remaining four radiogalaxies are in directions with lower exposure times. This correlation between nearby 
extended radiogalaxies and a subset of UHECRs is significant at the 99.9% level. For the remaining ~20 UHECRs, an isotropic 
distribution cannot be ruled out at high significance. The correlation found by the Auger Collaboration between the 27 UHECRs and 
AGNs in the Veron-Cetty & Veron catalog at D<71 Mpc has a much lower significance when one considers only the ~20 UHECRs 
not 'matched' to nearby extended radiogalaxies. No clear correlation is seen between UHECRs and supernovae, supernova remnants, 
nearby galaxies, or nearby groups and clusters of galaxies. 

Conclusions. Nearby extended radiogalaxies are the most likely source of at least some UHECRs detected by the Pierre Auger 
Observatory. The remaining UHECRs are not inconsistent with an isotropic distribution; their correlation to nearby AGNs is much 
less significant than earlier estimated. This is the first direct observational proof that radio galaxies are a significant source of UHECRs. 
The primary difference between the UHECR detections at the Pierre Auger Observatory and previous experiments, e.g. AGASA, may 
thus primarily be that the Southern Hemisphere is more privileged with respect to nearby radiogalaxies with highly extended radio 
jets and lobes. 
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1. Introduction 

Ultra-High Energy Cosmic Rays (UHECRs) are protons or light 
nuclei with energies greater than about 10 19 eV (10 EeV). When 
these highly energetic particles enter the earth's atmosphere 
they produce a shower of secondary particles and excite the 
Nitrogen molecules in the atmosphere; both effects can be de- 
tected from the ground. Recent advances in ground shower de- 
tectors allow a precise measurement of both the initial energy 
and the arrival directi on of the UH ECR above the earths atmo- 
sphere, e.g. AGASA jTakeda et ail 1 19991) and the Pierre Auge r 
Observatory dAbraham et all 12004 Ide Mello Neto et aUl2007l) . 
While UHECRs are rare they have a unique astronomical poten- 
tial: UHECRs with energies >50 EeV are not expected to be de- 
flected significantly by Galactic or intergalactic magnetic fields. 
Their measured arrival directions can therefore be traced directly 
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back to the originating source within the observatory measure- 
ment uncertainties. 

Suspected sources of high energy cosmic rays include core 
collapse Supernovae (SN II), Supernova Remnants, Pulsars, 
Active Galaxies - especially radio jets and radio lobes, and 
gamm a rays bursts ( for reviews and lectures see iKachelriesi 
(120081) ; iHillasI (I 19991) and references therein). In these sources, 
the cosmic rays are posited to undergo successive accelerations 
via scattering off energetic charged particles and/or in shocks. 
More exotic explanations include an origin in dark matter anhi- 
lation. 

The main energy loss for UHECR propagating over cos- 
mological distances is expected to be pion-pro duction, the so 
called Greisen Zatsepin K uzmin (GZK) effect dGreisenl [1966; 
IZatsepin & Kuzmiiull966l) . In this process, a UHECR interacts 
with a CMB photon and loses an estimated ~30% of its energy. 
Both the energy loss and the mean free path between interactions 
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depend on several details; energy l oss leng t h predictions vary 
betwee n 20 and 100 Mpc (see e.g. IStanevf d2004l)lKachelriessl 
(|2008j),RA08). If the GZK effect is at work then UHECRs with 
energy above ~50 EeV are expected to originate primarily in 
sources closer than this predicted energy loss length. 

The Pierre Auger Observatory (hereafter PAO) in Argentina 
dAbraham et al.L l2004t Ide Mello Neto et al.L 120071) . an array of 
1600 Cerenkov detectors spread over 3000 km 2 plus six op- 
tical telescopes, is designed to measure arrival directions and 
energies of cosmic rays via their secondary particle showers 
and their atmospheric flouresence. Latest details on the in- 
strumentation and exciting science results being produced by 
the Auger Observatory and Collaboration can be found at the 
Auger web-site http://www.auger.org The array has been in 
partial operation since January 2004. The Auger collabora- 
tion has reported 81 events between 2004 January 1 and 2007 
August 3 1 with reconstructed energies above 40 EeV and zenith 
angles smaller than 60°; of these, 27 have energi es above 
56 EeV dPierre Auger Collaboration et al.Ll2007ll2008l hereafter 
PA07,PA08). The origins of the latter 27 UHECRs are the focus 
of this paper. 

The Auger Collaboration found that the arrival directions of 
the 27 UHECRs with energies above 56 EeV are not isotropic 
at a 99% significance level (PA07,PA08). Instead, they found 
that the arrival directions are correla ted with the pos i tions of 
AGNs within -71 Mpc (PA07.PA08 lMollerach et ail d2008l) ). 
For this analysis they used all AGNs (and galaxies with HII 
nuclei) in the catalog of Ve ron-Cetty & Veron, 12th Edition 
dVeron-Cettv & Veron! |2006|). In their statistical correlation of 
all 81 UHECRs at energy above 40 EeV with the Veron-Cetty & 
Veron catalog of AGNs, PA07,PA08 found that the highest cor- 
relation between UHECRs and AGNs was obtained for a max- 
imum AGN distance of 71 Mpc, an angular separation of 3.2° 
between "matched" AGNs and UHECRs, and an energy thresh- 
old for UHECRs of 56 EeV. This maximum dista nce, 7 1 Mpc, is 
in line with the expectations of the GZK effect. Gorbunov et al.l 
(2007), in an astro-ph comment, disagreed with this finding, 
pointing out that weighting each AGN with its distance would 
predict a very different distribution of UHECRs than that found 
by the PAO. 

Our interest in this topic was provoked when a quick exam- 
ination of the 27 UHECR events showed that the direction of 
at least five of the 27 UHECR events were directly in the line of 
sight to nearby extended radiogalaxies. In this article, we attempt 
a more comprehensive correlation between the arrival directions 
of the 27 UHECRs with energies above 56 EeV and catalogs of 
potential sources of UHECRs, both Galactic and Extragalactic. 
Sec. [2] describes the data used in our study, Sec. [3] describes 
the principle results obtained, and Sec. [4] contains a brief dis- 
cussion and conclusions of our study. Distances to galaxies are 
calculated using 72 km s _1 Mpc -1 , except for relatively nearby 
galaxies for which we use distances as referenced. 

2. Data 

Positions and energies of the 27 UHECRs detected by the PAO 
with energies > 56 EeV were taken from PA08. The figures in 
this publication are slightly different from those of the Auger 
collaboration: PA07 and PA08 used the Aitoff projection (an 
equidistant projection with 2:1 aspect ratio), while we use the 
Aitoff-Hammer projection (an equal area projection). The lat- 
ter is more commonly used in astronomy and is typically re- 
ferred to as simply the 'Aitoff' projection, e.g. the IDL proce- 
dure aitoff.pro. Positions and energies of the UHECRs detected 



by AGASA are taken from Hayashida et al. (2000) which is an 
updation of the data listed in Take da et al.l ([T999). 

Our catalogs of astronomical sources are drawn from vari- 
ous publications and web-based source lists. In all cases we at- 
tempted to assemble as complete a sample as possible. Further, 
following the results of PA07, PA08, and the expectations from 
the GZK effect, we considered only objects within a distance of 
~200 Mpc. For matches to the PAO UHECRs we consider only 
those sources which rise above elevation 30° at PAO since the 
PAO UHECR list includes only UHECRs detected above this 
elevation. 

Our list of extragalactic supernova re mnants, galaxi e s, and 
galaxy clusters were drawn from NED (Ma doreetafl. |1992|) 
using the 'All Sky Search' option. The list of Gamma Ray 
Bursts (GRBs) was taken from the Gamma-ray Burst Real- 
Time Sky MapQ. For Galactic Supernova Remnants we used the 
list maintained by D. GreerQ- Extragalactic Radio Supernova 
were taken from Kurt Weiler's list at the web-site of the Naval 
Research LaboratorjQ. For extragalactic supernova, we used 
the Sternberg Astronomical Institute (SAI) Supernova Catalog 
maintained by Tsvetkov et al. In the case of radiogalaxies and 
radio s ources we used seve ral surveys and catalogs, incl uding 
NVSS dCondon et all 1 19981) . SUMMS dBocketall 1 19991) . and 
NED. The lis t of ex tragalactic jets is based on that compiled by 
Liu &Zhang| (120021) : to this list we added data on the galaxies as 
given in NED, and re-measured the total flux and total extent of 
the radio emission from NVSS and SUMMS maps. The total ex- 
tent of the radio emission therefore includes both radio jets and 
any radio lobes. This extent is typically referred to as the Largest 
Angular Size (LAS) or Largest Linear Size (LLS). In this work 
we added the LLS of the individual jets and lobes as scalars in- 
stead of vectors in order to discount the effects of jet and lobe 
bending. 

To recreate the results of the Auger Collaboration 
(PA07,PA08), we used the 12th edition of the 
IVeron-Cettv & Veronl (120061) catalog of AGNs. The list of 
AGNs (agn.dat) in this catalog includes nuclei with HII spectra 
which were previously misclassified as AGNs. These should 
ideally be left out of the correlation analysis. From the AGN 
statistics quoted in PA07,PA08, it appears that the Auger 
Collaboration included these HII nuclei in their analysis, so 
for consistency we do the same. Additionally, the list of QSOs 
(qso.dat) lists seven galaxies with redshift=0. Four of these are 
truly nearby galaxies, and the other three are probably high 
redshift quasars. Again, for consistency with PA07 and PA08 
we keep the latter three in our analysis. 

We have a lso used the Fanaroff-Ri ley (FR) classification for 
radiogalaxies dFanaroff&Rilevlll974 . This division is based on 
whether the jets increase (FR II) or decrease (FR I) in brightness 
with increasing distance from the nucleus. A clear division in 
radio luminosity is found between the two cl asses, with FR lis 
being more luminous in the radio than FR Is dOwen & Ledlowl 
1994). 



3. Results 

A comparison of the arrival directions of the 27 UHECRs and the 
p ositions of nearby galaxies with radio jets (all galaxies listed 
inlLiu & Zhangl |2002|) is shown in Fig.fl] We have divided the 
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galaxies with radio jets into three redshift bins: "nearby" (D < 
75 Mpc), "intermediate" (75 Mpc < D < 150 Mpc), and "dis- 
tant" (150 Mpc < D < 210 Mpc). In the closest redshift bin, D < 
75 Mpc, we distinguish between galaxies with radio structures 
more extended than 180 kpc ("extended") and those with radio 
structures less extended than 180 kpc ("compact"). We calcu- 
lated the total extent of the radio structure of any given galaxy 
by summing the lengths of the two radio jets and/or lobes as 
scalars instead of vectors. The total radio extent includes radio 
jets and any radio lobes. 

There is no intrinsic reason to use a radius ifi = 3.5° for a 
"match" between UHECR arrival direction and the position of 
an astronomical source. The instrumental resolution of the PAO 
is about 1° at these UHECR energies, to whic h one must a dd 
a small angle - perhaps up to a few degrees dStanevL 12004 - 
to account for the deflection of the UHECR by magnetic fields. 
The best correlation between UHECRs and nearby AGNs found 
by PA07, PA08 is based on considering if/ ~3.2°as a match. 
However, as we show later (e.g. Fig. [7]), the value of if/ in the 
analysis of PA07,PA08 should really be if/ ~ 3.1°^.3°. For the 
above reasons, we consider two values of if/ for a match: if/ = 1.5° 
and if/ = 3.5°. The distance limit for 'nearby' galaxies - 75 Mpc - 
is similar to the value d max ~ 71 Mpc found by PA07 and PA08. 
The division between 'compact' and 'extended' radio structures 

- 180 kpc - is the one 'new' free parameter in our analysis. We 
note that this value is close to the minimum linear size of the 
radio structures of a typical radio galaxy (e.i. lHardcastle et all 
fl998h . 

Of all 27 UHECRs with energy above 56 EeV detected by 
the PAO four (eight) are within if/ = 1.5° (if/ = 3.5°) of the ra- 
dio structures of nearby extended radiogalaxies in the field of 
view of the PAO. Conversely, of all ten nearby extended radio- 
galaxies in the field ov view of the PAO three (six) have ra- 
dio structures which can be matched to within 1.5° (3.5°) of a 
UHECR event. Of the remaining four 'unmatched' radiogalax- 
ies, three - NGC 1316, NGC 4261 and NGC 4760 - fall outside 
the area of maximum exposure at the PAO. For the intermediate 
and distant redshift bins there is little match between extragalac- 
tic radio jets and UHECRs; the only two matches within 3.5° are 
CGCG 403-019, a Seyfert/BL Lac at 1 12 Mpc, and Mrk 612, a 
Seyfert galaxy at 85 Mpc. 

The match between a UHECR and CGCG 403-019 de- 
serveres furter commen t. CGCG 403-019 is a BL Lac with 
a Seyfert 1 spectru m dVeron-Cettv & Veronl 119931) . In the 
12th AGN Catalog of IVeron-Cettv & Veronl d2006l) it is listed 
by its alias PKS 2201+04 and is classified as a Seyfert 1. 
It is a poweful radio source wi th a core- jet structure with 
total extent -6' or -190 kpc dUlvestad & Johnstonl Il984t 
lLaurent-Muehleisen et all 1 1993b . Within the unified scheme of 
radiogalaxies, BL Lacs are posited to be pole-on FR Is, implying 
a much larger deprojected radio extent. We therefore consider 
CGCG 403-019 as an 'intermediate' distance 'extended' radio- 
galaxy. Further, CGCG 403-019 is unique in being the nearest 
B L Lac with extended ra dio structure. The 12th AGN Catalog 
of I Veron-Cettv & Veronl d2006l) includes only three confirmed 
BL Lacs - V Zw 331, RXS J05055+0416, and TEX 0554+534 

- and only one probable BL Lac - RXS J2123 1-1036 - at D < 
150 Mpc. All four show compact radio emission in NVSS radio 
maps. 

The Auger Collaboration noted that the region around Cen A 
has an unusually large number of UHECRs; this is also the 
area with the maximum number of nearby extended radiogalax- 
ies. As seen in Fig. |2j the two UHECRs centered on Cen A 
can be explained as originating in either of Cen A, NGC 5090 



(D=48 Mpc; total radio extent -187 kpc), or even PKS 1308- 
441 (D=21 1 Mpc; total radio extent -1 140 kpc). The event just 
further south can be matched to the southern jet of Cen A, or to 
the relatively distant PKS 1308-441. CenB (D=54 Mpc; total ra- 
dio extent -250 kpc) is almost perfectly centered on a UHECR, 
and WKK 4552 (D=67 Mpc; total radio extent -360 kpc) is 
within -3° of another UHECR. The UHECR just to the north 
of Cen A is bracketed by IC 4296 (D=52 Mpc; total radio extent 
-1560 kpc) and ESO 443-G024 (D=71 Mpc; total radio extent 
-60 kpc). The very extended radio structure of the former actu- 
ally intersects with the 3.5° circle around the UHECR event. 

Are the nearby extended radiogalaxies matched to PAO 
UHECRs different from those not matched to PAO UHECRs? 
A comparison of the relevant properties of the ten nearest ex- 
tended radiogalaxies is given in Table 1, and their radio contour 
maps in Fig. [3] Table 1 also lists the relevant radio properties 
of CGCG 403-019 and Mrk 0612, the two intermediate distance 
AGNs which match to a UHECR. There is some difference in 
the combination of two factors: the radiogalaxies in the direc- 
tion to UHECRs include all of the subset which have radio jet 
morphologies closer to FR II, and are all in the area where the 
PAO has obtained a higher exposure time. In any case, given the 
small number statistics, i.e. typically 1 + 1 UHECR events per 
matched radiogalaxy, it is not surprising that UHECRs are not 
yet detected towards the remaining four radiogalaxies. Clearly 
more UHECRs are needed to form any firm conclusion. 

Fig. |4] compares the arrival directions of UHECRs with the 
distribution of nearby groups and clusters of galaxies. While 
there is an over-density of galaxies, clusters of galaxies and 
UHECRs around the region of Cen A, the overall distributions of 
nearby galaxies and nearby galaxy clusters do not correlate with 
the arrival directions of UHECRs. This argues against an ori- 
gin of UHECRs in material associated with lar ge mass haloes, 
e.g. D ark Matter. The more detailed analysis of iGorbunov et al.l 
(2007) also predicts a different distribution of UHECRs when 
the distances to galaxies and clusters are taken into account. 

We have also compared the arrival directions of the 27 
UHECRs with the positions of Galactic and nearby extragalac- 
tic supernovae and supernova remnants, of nearby extragalactic 
radio supernovae, of nearby BL Lacs, and of gamma-ray bursts. 
In all cases there appears to be no more than a random overlap. 

The analysis of PA07 and PA08 considered a single source 
population for the UHECRs. How would their results change if 
one were to consider multiple source populations, at least one 
of which is nearby extended radiogalaxies? Fig.|5]compares the 
histograms of the angular distribution of the 27 UHECRs with 
energies above 56 EeV, the subset of 19 UHECRs not towards 
nearby extended radiogalaxies, and the expectations from a ran- 
dom distribution of 27 UHECRs. The deviation from isotropy 
discussed by PA07 and PA08 (dashed line in the figure) comes 
primarily from the clustering of events around Cen A. Deleting 
the UHECRs which lie within 3.5° of nearby radiogalaxies (solid 
line in the figure) gives a distribution in which isotropy cannot 
be ruled out at high significance. 

The Auger Collaboration cross correlated the positions 
of the 27 UHECRs with th e AGNs in the 12th catalog of 
dVeron-Cettv & Veronl 2006) by varying the maximum AGN 
distance (z max ), the angular separation between AGN and 
UHECR to be considered a match (if/), and the threshold en- 
ergy of UHECRs. They then looked for the minimum cumula- 
tive probability of chance coincidences. Note that they appar- 
ently included HII nuclei from the catalog in their analysis. In 
this process (see PA08 for details) they found that the best match 
was found for z max = 0.017 (D - 71 Mpc) and if/ = 3.2°. We re- 
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peated their cumulative probability scans in the parameters z max 
and if/; we were unable to scan in energy as we do not have ac- 
cess to the position and energy data on the UHECRs with ener- 
gies less than 56 EeV. The resulting scans are shown in Fig. [6] 
The dashed line shows the result if all 27 UHECRs are used, i.e. 
the result obtained by PA08. In this case 19 of the 27 UHECRs 
match to some AGN or HII nuclei in the catalog. The middle 
solid line of each panel shows the result of considering only the 
23 UHECRs remaining after deleting the four UHECRs which 
best match (within 1.5°) the radio structures of nearby extended 
radiogalaxies (2 UHECRs for Cen A, and one each for Cen B 
and NGC 7626). In this case 13 of the 23 UHECRs match to 
some AGN or HII nuclei in the catalog. The upper solid lines 
show the results of considering only the 19 UHECRs remain- 
ing after deleting the four UHECRs mentioned above, plus the 
UHECRs matched within 3.5° to the South jet of Cen A, to 
CGCG 403-019, to WKK 4452, and to IC 4296. In this case 
1 1 of the 19 UHECRs match to some AGN or H II nuclei in the 
catalog. While the overall shape of the scans in if/ and z max re- 
main roughly the same, the significance of the results decrease 
significantly. In the case of the scan in if/, the minimum at around 
3.2° is no longer a clear and unique minimum. We repeated the 
above exercise after making two modifications to the AGN list: 
we deleted all galaxies with H II type nuclei (50 galaxies at z 
< 0.024) and the three probable high redshift QSOs which are 
listed as having z=0 in the AGN Catalog 12th. ed.. The results 
are shown in Fig. [7] The main difference with the previous figure 
is the broader valley in the if/ scan. Instead of a clearly defined 
minimum at if/ — 3.2°one now has a greater range if/ ~ 3.2°- 4.3°. 

4. Discussion & Concluding Remarks 

The arrival directions a subset of the 27 UHECRs with ener- 
gies above 56 EeV detected by the PAO are statistically most 
closely related to nearby radiogalaxies with extended radio jets 
and/or lobes. It is thus likely that nearby (due to the GZK ef- 
fect) radiogalaxies with extended (radio extent > 180 kpc), two 
sided jets and lobes are the source population of a subset of the 
UHECRs detected by the PAO. Additionally, there is weak evi- 
dence that hosting a jet with morphology closer to FR II makes a 
radiogalaxy more likely to be a UHECR source. Interestingly, all 
these factors are consistent with the lack of UHECR detections 
towards M 87. 

The results of the previous section were based only on 
UHECRs detected by the PAO. We have also tested for corre- 
lations against the UHECRs with energy above 56 EeV detected 
by AGASA (red circles with radius 3.5° in FigQ] Among the 
nearby (D < 75 Mpc) galaxies with extended radio jets, the only 
clear match is to CGCG 514-050, a radiogalaxy at D=72.2 Mpc 
with radio extent ~250 kpc. The remaining three nearby galax- 
ies with extended radio galaxies, NGC 315, NGC 383, and 
NGC 5127, are not matched to within 3.5° of any AGASA 
UHECR event, though there are several UHECRs at slightly 
larger distances. Additionally, One AGASA UHECR is close to 
the location of NGC 7626 and another is centered on the radio- 
galaxy 3C 120, at D=135 Mpc and total radio extent ~3 15 kpc. 

Among the UHECRs detected by PAO, an isotropic distri- 
bution of the ~20 UHECRs which are not in the line of sight 
to nearby extended radiogalaxies cannot be ruled out at high 
significance. The significance of the correlation found by PA07 
and PA08 between UHECR arrival directions and nearby AGNs 
from the catalog of Veron-Cetty & Veron is much lower when 
one only considers the ~ 20 UHECRs not matched to nearby 
extended radiogalaxies and even lower when H II nuclei in the 



the AGN Catalog are deleted. The findings of PA07,PA08 were 
probably largely influenced by the correlation between UHECRs 
and nearby radiogalaxies which we have shown above. 

The UHECR arrival directions are not strongly corre- 
lated with either supernovae, extragalactic radio-supernovae, or 
nearby groups and clusters of galaxies. 

The main difference between the results of the Pierre Auger 
Collaboration and those from previous studies of UHECRs, e.g. 
AGASA, could primarily be that the southern location of the 
PAO is more privileged with respect to nearby extended radio- 
galaxies. 
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Fig. 1. A comparison of the arrival directions of the 27 UHECRs with the positions of galaxies hosting radio jets. PAO UHECR 
events with energies above 56 EeV are shown as black circles with radius 3.5°. AGASA UHECR events with energies above 56 EeV 
are shown as red circles with radius 3.5°. The colour symbols mark the positions of galaxies with radio jets at D < 75 Mpc (red), 
75 Mpc < D < 150 Mpc (blue), and 150 Mpc < D < 210 Mpc (black). For galaxies at D < 75 Mpc, filled symbols are used for 
galaxies with radio structures more extended than 180 kpc, and open symbols for galaxies with radio structures less extended than 
180 kpc. Galaxy names are marked in some of the less crowded regions in the plot, in the same color as the corresponding symbol. 
The SuperGalactic plane is marked by the dashed line. In this and following figures, we use an Aitoff Hammer Projection. The 
orientation and labelling of the longitude axis follows PA07 and PA08 for easier comparison. 
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Fig. 2. A detail of Fig. Q] showing the region around the radiogalaxy Cen A. Symbols are the same as in Fig. Q]except in the case of 
Cen A, for which we use a cross to mark the galaxy center and a polygon to outline its radio jet structure as seen in a 403 MHz map. 
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Fig. 3. The radio morphologies of the radiogalaxies in Table 1. Unless otherwise specified, the radio maps are from NVSS. From 
left to right, the contour maps are top row: NGC 1316 (SUMMS), NGC 3557, NGC 4261, NGC 4760; middle row: NGC 5090 
(SUMMS), Cen A (408 MHz), IC 4296, Cen B (SUMMS); bottom row: WKK 4452 (SUMMS), NGC 7626. Radiogalaxies which 
match the positions of UHECRs (middle and bottom rows) are preferentially those with intermediate FRI-FRII type characteristics. 
Those not matched to positions of UHECRs (top row) are FR Is or in the case of NGC 1316, a 'relaxed double'. 
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Fig. 4. A comparison of the arrival directions of the 27 UHECRs with the positions of clusters of galaxies (as listed and defined 
by NED). The symbols for UHECRs are the same as in Fig.Q] Red symbols are used for galaxy clusters at D < 75 Mpc, and blue 
symbols for galaxy clusters at 75 Mpc < D < 200 Mpc. 
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Fig. 5. Histograms of the angular separation between UHECRs detected by the PAO. The dashed line is for all 27 UHECRs with 
energy above 56 EeV. The solid line is for the subset of 19 UHECRs which are not matched to nearby radiogalaxies with extended 
radio structures (see text). The dotted line is the expectation for a uniform distribution of the 27 UHECRs given the exposure time 
weighting in declination for the PAO. The histograms are offset by a few degrees in x for clarity. 
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Fig. 6. Probability of the null hypothesis (isotropic distribution) vs. maximum angular distribution \p (left panel) and maximum 
AGN redshift z max (right panel). In thi s figure we followed th e analy sis of PA07,PA08, i.e. correlating UHECRs to all galaxies at z 
< 0.024 in the AGN Catalog 12th ed. dV6ron-Cetty & Veronl [20061) . The dashed lines are the results of using all 27 UHECRs, i.e. 
the same plots as in PA07 and PA08. The middle solid line is the result of deleting all 4 UHECRs within 1.5° of the radio structures 
of nearby extended radiogalaxies and the upper solid line is the result of deleting all 8 UHECRs within 3.5° of the radio structures 
of nearby extended radiogalaxies (see text). 
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Fig. 7. The same as Fig.|6j but this time, of all galaxies listed in the AGN Catalog 12th ed., we do not consider galaxies with H II 
type nuclei and the three high redshift QSOs which are listed as z=0 sources (see text). The main difference from Fig.|6]is that the 
minimum value for the null hypotesis in tfr is now over a broader range: tfr ~ 3.1°to 4.3°. 
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Table 1 



Object 


1 


b 


Dist. 


S4.85GHZ 


S 1.4GHz 


Ss43MHz 


S408MHZ 


Size 


Morph. 


Ref. 


name 


(Degrees) 


(Degrees) 


(Mpc.) 


(Jy) 


(Jy) 


(Jy) 


Jy 


(kpc.) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC 1316 (Fornax A) 


240.162919 


-56.690021 


24.5 






169 


259 


375 


Rel.Dbl 


2,3 


MRK 0612 


187.717486 


-45.039107 


85.2 




0.013 






2.6 


Sey.2 


8 


NGC 3557 


281.578358 


21.089037 


42.5 


0.32 


0.63 




0.96 


310 


FRI 


2 


NGC 4261 (3C 270) 


281.804700 


67.372841 


31.1 


4.6 


19.4 




41.5 


436 


FRI 


2 


NGC 4760 


303.609402 


52.375264 


66.2 


0.37 


0.9 




4.7 


322 


FRI 


5,7 


NGC 5090 


308.611065 


18.838673 


47.6 


1.98 


5.8 


7.16 




187 


FRI 


3 


NGC 5128 (CenA) 


309.515874 


19.417325 


3.4 


62.84 


278 


392 


2746 


453 


Int. FRI/II 


3 


IC4296 


313.538349 


27.973003 


51.9 


1.88 




14.9 


34 


1564 


Int FRI/II 


2,3 


PKS 1343-60 (CenB) 


309.721407 


1.731590 


53.8 


27 


79 






256 


FRI 


6 


WKK4452 


321.594046 


6.291192 


66.7 




3.2 


3.8 




358 


FRI+FRU 


1,3 


CGCG 403-019 


64.701023 


-38.631307 


112.5 


0.74 


0.78 




1.64 


187 


BLLac/Sey.l 


4 


NGC 7626 


87.859111 


-48.378945 


47 


0.25 


0.86 




1.9 


212 


FRI 


2 



Note. — This table lists some relevant properties of all ten nearby and extended radiogalaxies (see text) in the 'field of view' (see 
text) of the Pierre Auger Observatory, plus two Active Galaxies which potentially match to UHECRs detected by the Pierre Auger 
Observatory. 

Note. — Columns are: (1) galaxy name; (2) and (3) Galactic longitude and latitude; (4) distance; (5) to (8) approximate total flux in low 
resolution radio maps, taken from NED, measured from NVSS or SUMMS, or from the references in Col. 1 1; (9) largest linear size, in kpc, of 
the radio structures; (10) jet morphology (see references in Col. 11). Here FRI/II is intermediate between FR I and FR II, FRI+FRII implies 
one jet has an FR I morphology and the other a FR II morphology, Rel. Dbl. is a 'relaxed double', and Sey. denotes Seyfert. (11) References. 
(1) Jones 1986; (2) Birkinshaw & Davies (1985); (3) Jones & McAdam; (4) Ulvestad & Johnston (1984); (5) Saikia et al. (1986); (6) Jones et 
al. (2001); (7) NVSS; (8) Nagar et al. (1999). 



